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SYSTEM AND PROCESS FOR MEASURING, COMPENSATING AND TESTING 
NUMERICALLY CONTROLLED MACHINE TOOL HEADS AND/OR TABLES 

CROSS-REFERENCE TO RELATED APPLICATION 

5 The present application is a United States National Phase Application of International 

Application PCT/IT2003/000556, titled "System and Process for Measuring, Compensating 
and Testing Numerically Controlled Machine Tool Heads and/or Tables," filed September 19, 
2003, which claims priority from International Application PCT/IT02/0065 1 , titled "System 
and Process for Measuring, Compensating and Testing Numerically Controlled Machine Tool 
10 Heads and/or Tables," filed October 11, 2002, the contents of which are incorporated herein 
in their entirety. 

The present invention relates to processes for measuring, compensating and testing 
rotary heads and/or tables for numerically controlled machine tools and to a hardware and 
software system that performs such processes. 

15 The invention will be described and shown herein below as applied to machine heads, 

but obviously everything that is described can be similarly applied to systems with head + 
table or only machine tool tables obtaining the same inventive effects and advantages, since 
simple Cartesian reference systems conversions are hereby involved. 

The rotary heads (herein below called "heads") of numerically controlled machine tools 

20 are machine tools axes controlled through a numeric control (CNC). These heads are 

commonly equipped with one or two rotation degrees of freedom. Possible variations to such 
configuration, that has to be deemed included within the scope of the present invention, are 
rotary heads in which the rotation movements are partially or wholly decomposed on rotating 
tables. 

25 The hardware system, that will be described herein below, allows measuring with 

automatic processes enough values so that measuring processes and software can process 
compensations and tests of this part of the machine. 

The term "geometric measure" of the heads means the measure of geometric head 
errors. These errors are generated by an incorrect component assembling, by an incorrect 
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realization of the components or by the component weights. Herein below a more detailed 
classification of these error types will also be explained. 

The term "compensation" of the heads means the software correction of the measured 
errors within the numeric control. As previously mentioned, an error classification will be 
5 given below and therefrom the related software compensation algorithms will follow. The 
compensation further includes both a measuring software and a compensation software. The 
measuring software allows performing and processing the measuring data while the 
compensation software performs corrections of measured and processed errors. 

The term "automatic processes" of the heads means a series of automatic machine 

10 movements that, through the hardware system, collect measures that the measuring software 
processes to create compensation software input data. 

The term "dynamic test" of the heads means the measure of dynamic responses of 
machine movements that involve rotation head axes. Typically, these tests point out 
"twitching" and "ripple" problems for these axes. 

15 Currently, there are no dynamic test systems for the heads, while there are manual or 

semi-automatic geometric measuring processes: these latter ones, above all for the purposes of 
the present invention, are not systems that are integrated in the numeric control and they do 
not allow the freedom to suit at will the complexity of the error model that has to be described 
by performing the measures. 

20 The current geometric measures provide for the use of measuring systems, such as 

decimal or millesimal comparators, and artefacts, such as certified squares, control cylinders 
(commonly called "gage tools"), etc. The measure operator, with a series of known manual 
operations, after having assembled and repeatedly positioned the above mentioned artefacts, 
takes note of the measure values. The semi-automatic processes mainly use digital feeler 

25 devices (see, for example, the Renishaw system MP 10) and a reference artefact (mainly a 
ball). With measuring processes it is thereby possible to obtain part of the values that can be 
measured with the manual system. The major problem with these known techniques is that 
they make the integration with the numeric control difficult and make the use of complex 
compensation models unfeasible, since the number of measures to be performed and the timing 

30 generated by such technique would be excessively long and cumbersome. 
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One of the limits of such applications is that not all the necessary measures and/or 
variables for the geometric head measure can be automatically measured, as can be well seen 
from the schematic drawings in Fig. 1. In such Figure, it can be seen that the parallelism 
check between the axis of head 1 equipped with two parts 3 and 5 and a stem 7 in a mutual 
5 rotation movement, and the reference Cartesian axes is carried out through a measuring 

instrument (commonly called "gage tool") 9 placed next to the stem 7. The limits of such type 
of operation essentially are the adopted measuring systems and the following limited measuring 
processes. 

The object of the present invention is solving the prior art deficiencies, by providing a 
10 system and a process for measuring, compensating and testing numerically controlled machine 
tool heads that allow obtaining the following innovations: 

a) an instrument, and therefore processes, for automatically performing the geometric 
measures integrating such system into the numeric control; 

b) the chance of increasing the complexity of such measures in order to describe more 
15 complex head errors, always having, however, measure execution times that are extremely 

reduced with respect to the prior art; 

c) the use of such instrument for dynamic tests; 

d) the use of known correction systems to perform the compensation of measured 

errors. 

20 The above and other objects and advantages of the invention, as will appear from the 

following description, are reached by systems and processes as disclosed in the independent 
Claims. Preferred embodiments and non-trivial variations of the present invention are 
disclosed in the dependent Claims. 

The present invention will be better described by some preferred embodiments, 
25 provided as a non-limiting example, with reference to the enclosed drawings, in which: 

Figure 1 is a schematic view that shows the measuring steps of axis quadrature in the 
prior art; 

Figures 2A to 2D are views of the support base 1 1 part of an embodiment of the system 
of the present invention; 
30 Figure 3 is a side view of the gage tool part of an embodiment of the system of the 

present invention; 



-3- 



Figures 4A to 4C are schematic views of an embodiment of the head of the system of 
the invention with related geometric models; 

Figures 5A to 5C are schematic views of another embodiment of the head of the system 
of the invention with related geometric models; 
5 Figure 6 is a schematic view of an embodiment of the system of the invention; 

Figures 7 A to 7C are schematic views that exemplify the way in which the measured 
amounts are shown; 

Figures 8A to 8C are schematic views that show the movements of the support base 1 1 
of Fig. 2A to 2D; 

10 Figure 9 is a schematic view that shows an application example of the system of the 

invention; 

Figures 10 to 13 are schematic views of further steps in the example in Fig. 9; 

Figures 14A and 14B are schematic views that show a further application example of 
the system of the invention; 
15 Figure 15 is a perspective view that shows some operating positions of the system of 

the invention; 

Figure 16 is a schematic view that shows the nutator head; 

Figures 17 to 21 are schematic views of the movements of the system of the present 
invention; and 

20 Figures 22 to 24 are vector graphs that show the actions of the head model of the 

present invention. 

With reference to the Figures, a preferred, but non-limiting embodiment of the system 
of the present invention is shown. The system will be described herein below as applied to the 
field of measuring, compensating and testing numerically controlled machine tool heads 1, but 

25 it is obvious that it can find a valid and efficient application to any field in which a completely 
automated, accurate and reliable measure of object positions in reference systems is necessary. 
With reference to the Figures, and in particular to Fig. 2 to 5, the system for measuring, 
compensating and testing numerically controlled machine tool heads 1 and/or tables 
substantially comprises: 

30 at least one support base 11 equipped with a plurality of distance sensors 14; 

at least one device 16 of the gage tool type composed of an elongated cylinder 17; the 
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cylinder 17 is equipped at one of its ends with connection means 18 for the heads 1 and is 
equipped at another opposite end with a ball 20, that is placed next to the sensors 14 so that 
they are able, always and in any position, to measure the distance that separates them from the 
ball 20. 

5 In particular, the support base 1 1 can be of a circular shape and is preferably equipped 

with three distance sensors 14 placed on the base in positions that are mutually offset by 120°. 
Instead, the connection means 18 are of the tapered type and the heads 1 are adapted to 
receive, in one of their moving parts 3, 5, the connection means 18 for the unmovable 
connection thereto during the measures. 

10 According to what is shown in Fig. 6, the system of the invention is operatively 

coupled to processing means 30 comprising means 31 for performing processes for measuring 
errors that can be modeled, means 32 for performing processes for measuring errors that 
cannot be modeled and means 33 for performing dynamic checks. Moreover, the heads 1 are 
operatively coupled with test and control means 34 also comprising means 35 for performing 

15 compensation processes of errors that can be modeled and means 36 for performing 
compensation processes of errors that cannot be modeled. 

According to the measure of positioning errors of the ball 20 of the gage tool 16, that 
ideally shows the tool tip, two possible processes can be obtained: 

1) The head 1 is moved by using already active compensations (those of errors that can 

20 be modeled). Then the CNC deems to have kept the ball 20 centre unmoved (in CNC systems 
this type of movement is commonly called movement with enabled RTCP). The error measure 
can be carried out using one of the following alternatives: 

a) The values provided by the sensors 14 of the support base 1 1 can be transformed 
into real coordinates of the ball 20 of the gage tool 16. The difference between theoretical 

25 coordinates and real coordinates of the ball 20 shows the measured errors. 

b) The values provided by the sensors 14 of the support base 1 1 can be used to 
correct the position of the machine linear axis in order to take back the ball 20 to the point 
where the sensors 14 provide the initial values. In this way the ball 20 will not have been 
moved but otherwise the linear axis will have performed an additional movement with respect 

30 to the one that the CNC would have imposed to them depending on currently active 
compensations. Such corrections are the measured errors. 
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2) The head is handled without using any active compensation in the CNC. Through 
the previously-listed process b, the linear axes position is corrected in order to take back the 
ball 20 to the point in which the sensors 14 provide the initial values. Such induced movements 
represent the indirect error measure for computation purposes. 
5 As regards the error detection mechanism, through a software that can be easily and 

readily realized that performs some transforming, the support base 11 with the sensors 14 
provides the position of the ball 20 of the gage tool 16. The process to detect the position of 
the ball 20 is known in the art and will not be described in further detail in this context. For 
measuring the errors, it is necessary to perform a calibration process of the sensors 14 of the 

10 support base 1 1 (such process is known in the art and will not be described herein in detail). 
Independently from the technique used to extract the positioning errors from the 
controlled movements to the machine, as shown in Fig. 7 A to 7B, the measuring system is 
idealized assuming that the machine displacement (desired displacement + measured errors) is 
able to be measured. In Fig. 7A, a displacement error A of the head 1 axis is shown, 

15 depending on which, as shown in Fig. 7B, a measured error D occurs following the 

programmed movement of the head 1 from a position B to a position C, such movement being 
represented by the corresponding arrows in the Figure. The situation in Fig. 7C thereby 
occurs, that shows the diagram of the performed measure. 

As regards the error measuring process, the advantages of the configuration of the 

20 system with sensors 14 and support base 1 1 as a tripod are: 

1) Enough measuring accuracy in order to determine the centre of the ball 20 of the 
gage tool 16 and therefore the errors to be measured. The sensors 14 can be realized using 
different technologies and can provide or not the contact with the ball 20 of the gage tool 16. 
When distance sensors 14 are used instead of the contact ones and of a capacitive type, the 

25 support base 1 1 can be shaped as a hollow spherical cover (not shown) whose focal point is the 
centre of the ball 20 of the gage tool 16. This is in order to minimize capacitive measures 
noises. Such modifications (and similar ones) must not be deemed as variations from the basic 
concept of the support base 11 as previously described. 

2) Possibility of slanting the gage tool 16 below the horizontal line without collisions 
30 between moving bodies and measuring instrument. The wide range of positions that can be 
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reached for the measures allows the measuring software to inspect thereby all interesting head 
positions. 

According to what is shown in Fig. 8 A to 8C, the support base 1 1 can be connected to 
bearing means 40 adapted to allow a rotation of the support base 1 1 itself up to 90° with 
5 respect to its own axis (F), in order to reach a plurality of operating positions between two 
mutually perpendicular extreme axes (F, G). The bearing means 40 are further adapted to 
simultaneously allow a rotation of the support base 1 1 , once having reached the extreme axis 
(G) position, around the axis (F) perpendicular thereto. 

As regards the geometric measure process of the head 1, as described previously, the 
10 automatic measures allowed by the measuring instrument and the measure processes, allow 
increasing the complexity degree of the geometric head model. In this context, the term 
"geometric model" of the head 1 means the mathematical model that describes the real 
behaviour of the head 1 with respect to the theoretical one. The parameters of such model are 
those that are obtained through the measures. 
15 Fig. 4 A to 4C and 5 A to 5C show some types of heads and the related simplified 

geometric models that are used today related to a possible more complex model that can be 
used with the system described here. 

The parameters of these geometric models can be called errors that can be modeled, 
meaning that a possible compensation system is able to use these parameters to compute the 
20 errors that the CNC must correct for working purposes. These latter errors are, for example, 
positioning errors of the tool end (tool tip). 

When the head 1 has been compensated by using the above described system, it will 
have a series of errors that the used model, though more complex, does not describe. These 
errors, that can be called errors that cannot be modeled, result into a still incorrect tool tip 
25 positioning. The instrument described here (but also equivalent instruments such as digital 
feeler pins and reference balls) is also able to perform the measure of these residual errors. 

It will now be demonstrated that the described system is able to provide necessary data 
for identifying the parameters of a geometric model of a head 1 selected as an example. The 
demonstration is divided into three examples that describe geometric models with increasing 
30 complexity for measuring the errors that can be modeled and a further example that describes 
the measure of the errors that cannot be modeled. 
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For better clarifying, the examples of procedures will include: 
Example 1) Procedure for measuring a head with a simple geometric model 
Example 2) Procedure for measuring a head with a complex geometric model 
Example 3) Procedure for measuring a head using an algorithm for solving complex 

5 models 

Example 4) Procedure for measuring errors that cannot be modeled 
Example 5) Description of a practical arrangement 
EXAMPLE 1: 

Let us take into account the geometric model described in Fig. 9 where it is assumed as 
10 an example that the axes runs are A [0;360] B [-90;90] (refer to Fig. 9): 
The parameters of this model are: 

• Angular positioning accuracy for axes A and B. 

• Direction and rotation centre of axes A and B. The exact orientation of the rotation 
axis must be located in order to determine the parallelism position between rotation 

15 axis and linear axis. 

• Zero pre-set of axes A and B. The 0 pre-set shows the position in which an axis is 
programmed in position 0. For a head the 0 pre-set is the point in which the tool is 
aligned with axis Z. 

• P (commonly called pivot), DY (distance between rotation axis for axis B and axis 
20 A), DTy (distance in plane YZ between gage tool axis and rotation axis for axis A), 

DTx (distance in plane XZ between gage tool axis and rotation axis for axis A) 
It must be noted that the figure in plane YZ assumes that the 0 pre-set for axis A has 
already been located. Otherwise the distance between the rotation axes would have, in addition 
to a component DY, also a component DX since the rotation axis of axis B would not be 
25 perpendicular to plane YZ. It is for such reason that before computing DY, the 0 pre-sets are 
computed. 

For easy of exposure in the figure the rotation axes are shown parallel to the linear 
axis. In practice this is obtained after having measured and corrected the parallelism between 
rotation axes and linear axis. 
30 The measuring procedure consists in the following steps: 

Step 1: 
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Head positioning A=0 B=90°. 

Movement of axis A from 0° to 360° and positioning errors acquisition for a defined 
pitch in axis A movement. 

A series of points is then measured that describe a circle. The circle equation can be 
5 computed with known quadratic error minimizing methods. 

From the measured points (which have an univocal relationship with axis A transduced 
position) it is possible to determine the angular positioning accuracy for axis A. 

From the circle equation it is possible to determine the rotation axis slant for axis A 
with respect to axis Z (parallelism between rotation axis A and linear axis Z). 
10 From the measured points it is then possible to determine which is the axis A 

transducer value for which the measured points lie in plane YZ. The 0 pre-set for axis A is 
thereby determined. 

It must be noted that the errors in axis B (0 pre-set, parallelisms and transducer 
linearity) have no influence whatsoever. 
15 Step 2: 

With a similar procedure to the one in step 1, the head is positioned A=0 B = 90° and a 
movement of axis B is carried out from 90° to -90° and a positioning error acquisition is 
carried out for a defined pitch in axis B movement (see Fig. 15). 

A series of points is then measured that describe a semicircle. The circle equation can 
20 be computed with the known quadratic error minimizing methods. 

From the measured points (which have an univocal relationship with axis B transduced 
position) it is possible to determine the angular positioning accuracy for axis B. 

From the circle equation it is possible to determine the rotation axis slant for axis B 
with respect to axis X (parallelism between rotation axis B and linear axis X). 
25 From the measured points it is then possible to determine which is the axis B transducer 

value for which the measured points lie in plane XZ. The 0 pre-set for axis B is thereby 
determined. 

At this time both for axis A and for axis B the following are known: angular 
positioning accuracy, parallelism between rotation axes and linear axis, zero pre-set. By 
30 activating such parameters in the compensation software or by taking into account their value 



-9- 



in the measuring software, it can be assumed that their effects have been cancelled and thereby 
the head has been driven to the condition described in the starting figure. 
Step 3: 

Errors in positions A=0 B=0 and A = 180 B=0 are measured (see Fig. 10). 
The relationship: 2*DTx = X is obtained, 

where X is the index of a value that can be measured through the gage tool (see Fig. 3) 
+ tripod (see Fig. 2) measuring system. 
Step 4: 

Errors in positions A=0B=0 and A = 180 B=0 are measured (see Fig. 11). 
The relationship: 2*(DY+DTy) = Y is obtained, 

where Y is the index of a value that can be measured through the gage tool + tripod 
measuring system. 
Step 5: 

Errors in positions A=0 B = -90 and A ==180 B = 90 are measured (see Fig. 12). 
The relationship: 2*DTy = Z is obtained, 

where Z is the index of a value that can be measured through the gage tool + tripod 
measuring system. 
Step 6: 

Errors in positions A=0 B=-90 and A = 180 B = -90 are measured (see Fig. 13). 
The relationship: 2*(DY + P) = Y is obtained, 

where Y is the index of a value that can be measured through the gage tool + tripod 
measuring system. 

The equations found in steps 3 to 6 are all independent and their resolution is 
mathematically able to be demonstrated. 

With a similar procedure a nutator head can be measured. 

Fig. 16 shows a simplified geometric model of a generic head 1 (nutator head). It is 
assumed as an example that the axes runs are A [0;360] B [-180;180] 
The parameters (or errors to be measured) in this model are: 

• Angular positioning accuracy for axes A and B. 

• Direction and rotation centre for axes A and B. 

• Zero pre-set for axes A and B. 
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• DTx, DX, DTy, DY and P. The angle (a) of 45° is assumed as known. The 
measures to be carried out are 5 and the support base 1 1 settings are always those 
of the basic configuration. 
As in the previous case in Fig. 16, the rotation axes are shown as parallel to the 
5 theoretical configuration (axis A parallel to axis Z and axis B in plane XZ at 45° from axis X). 
In practice this is obtained after having measured and corrected the parallelism between 
rotation axes and linear axis and obtained the 0 pre-sets. 

The measuring procedure consists in the following steps: 
Step 1: 

10 As in the previous case, a circle is performed for measuring axis A. B is placed at 180° 

(gage tool in plane XY) and the measure is performed by moving A from 0 to 360 degrees. 

The operations to be performed are the same ones and it is thereby obtained: angular 
positioning accuracy for axis A; parallelism between rotation axis A and axis Z; 0 pre-set for 
axis A for which the gage tool lies in plane XZ. 

15 Step 2: 

As in the previous case, a circle is performed for measuring axis B. A is placed at 0° 
and the measure is performed by moving B from -180 to +180 degrees. 

The series of points describe a circle that lies on a plane that ideally should intersect 
plane YZ in a straight line parallel to Y and plane XZ in a straight line slanted by -45° with 
20 respect to X. Deviations from such configuration will be parallelisms corrections for axis B. 

As previously the corrections are detected in order to correct the angular positioning 
accuracy for B and the 0 pre-set for which the gage tool is vertical. 

Like in the previous example, being known both for axis A and for axis B transducer 
linearity, parallelism between rotation axes and linear axis and zero pre-set, the head has been 
25 driven to the condition described in Fig. 16. 

Step 3: 

In Fig. 17, the following is applied: 
Movement: 
AO = > 180 
30 B0 

Relationship: 
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DTx + DX = X 
Step 4: 

In Fig. 18, instead, the following is applied: 
Movement: 
5 AO 

BO => 180 
Relationship: 
P - DTx = X 
Step 5: 

10 In Fig. 19, the following is applied: 

Movement: 
AO = > 180 
B180 

Relationship: 
15 2P + 2DX = X 
Step 6: 

In Fig. 20, the following is applied: 
Movement: 
AO = > 180 
20 BO 

Relationship: 
2DY + 2DTy = X 
Step 7: 

In Fig. 21, the following is applied: 
25 Movement: 
AO 

BO => 180 
Relationship: 
2DTy = Y 
30 P = Z 
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The equations found in steps 3 to 7 are all independent and their resolution is 



mathematically able to be demonstrated. 

If angle a shown in the Figures is not 45° and is not known with accuracy, by using 
the analytical method described below, such angle would be measurable. 



Figures 22 to 24 show an advanced geometrical model of a generic head 1 . The 
parameters (or errors to be measured) of this model are: ax, Px, DBx, DTx, P, at, ay, py, 
DBy, DTy, pt. 

In these models the non-parallelisms are pointed out for the rotation axes with respect 
10 to linear axis (ax, px, ay, py). Such parameters are measured with the same procedures 

mentioned in steps 1 and 2 which further allow obtaining the transducer linearity and the 0 pre- 
set value. 

Having measured such parameters and by activating the compensation software or 
taking into account such values within the measuring software, it can be ideally assumed to 



15 have corrected such errors and have transformed the affected head into a simpler model (see 
Fig. 14A). 



20 thereby performing steps 3 to 6 plus other two steps that describe independent positions and 
equations, again a series of equations is found that are arithmetically able to be solved. It must 
be noted that for this model DTx + DX is a single variable (from now on DX). 



5 



EXAMPLE 2: 



Or, according to which is the reference model, see the example in Fig. 14B. 
As an example refer to the case in the second figure. 

With respect to the model in example 1 , we have parameters at and pt more and 



25 



30 



The steps therefore are: 

Step 3: positions A=0 B=0 and A=180 B=0 

Relationship: 2*(DX + P*sin(at)) = X 

Step 4: positions A=0 B=0 and A = 180 B=0 

Relationship: 2*(DY + DTy + P*sin(pt)) = Y 

Step 5: positions A=0 B = -90 and A = 180 B = 90 

Relationship: 2*(DTy + P*sin(pt)) = Z 

Step 6: positions A =0 B = -90 and A = 180 B = -90 

Relationship: 2*(DY + P*cos(P0) = Y 
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Step 7: positions A=0 B = -90 and A=0 B=90 

Relationship: 2*(P*sin(pt)) = Y 

Step 8: positions A=0 B=0 and A=0 B=-90 

Relationship: P*sin(Pt) - DTy - P*cos(pt) = Y 

P*sin(pt) + DTy + P*cos(pt) = Z 

Summarizing the processes disclosed in Examples 1 and 2, it is still possible to locate 
the following macro-steps for a generic head 1 : 

1. for a head (20), there are planes in which the circumference described by the tool tip 
lies when the following movements are realized: 

- axis B at 90°; axis A that performs one revolution (circle 1) 

- axis A at 0°; axis B that performs one revolution (circle 2), 

- these planes are perpendicular and parallel to machine tool Cartesian axes; such 
process comprises the steps of: 

- executing the circle 1 ; 

- rebuilding with mean square methods or the like the circle 1 through a series of 
points describing it; 

- locating the non-parallelism of the plane passing through the circle 1 with respect to 
the plane that is orthogonal to the ideal rotation axis of axis A; 

- locating the relationship between position transduced by axis A measuring systems 
and related tool tip location point and then computing the angular positioning 
accuracy of axis A; 

- executing the circle 2; 

- rebuilding with means square methods or the like the circle 2 through a series of 
points describing it; 

- locating the non-parallelism of the plane passing through the circle 2 with respect to 
the plane that is orthogonal to the ideal rotation axis for axis B; 

- locating the axis A position 0 through the component of the previously described 
angles that lies in the plane perpendicular to the ideal rotation axis of axis A; 

- locating the axis B position 0 through the measured point on circle 2 that allows 
having the tool as vertical; and 
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- locating the relationship between position transduced by axis B measuring systems 
and related tool tip positioning point and then computing the angular positioning 
accuracy of axis B 

2. According to the geometric model of the head 1 that has to be described, and 
5 therefore to the number of parameters that are still not known, a sequence of positionings of 
the head 1 is performed, adapted to obtain algebraic equations that link the above parameters 
to the measured errors; such equations must be linearly independent and must be equal to the 
number of parameters to be determined. 

EXAMPLE 3: 

10 The geometric model complexity could further increase (therefore increase the number 

of model parameters) till the search of solutions through algebraic equations that describe 
independent positions is made difficult. In this case more complex mathematics and more 
sophisticated solution algorithms could simply be adopted. 

More complex models are, for example, the Rodriguez-Hamilton models that describe 

15 the 6 degrees of freedom for coupling two bodies. In this case our error model would have 6 
degrees of errors for axis A movement, 6 for B axis movement and 6 for gage tool rotation in 
the spindle. 

In order to solve such type of problem, known techniques can be adopted that solve 
mathematical problems or models depending on a series of experimental measures (our error 

20 measures in the different head positions). These techniques are based on minimizing mean 

square errors or on using neural algorithms (as known in the art, for example, from: Kim K., 
Kim M.K. "Volumetric Accuracy based on Generalized Geometric Error Model in Multi-Axis 
Machine Tool", Mec. Mach. Theory, Vol. 26 (1991) No. 2, pages 207-219; Duffie N.A., 
Yang S.M., "Generation of Parametric Kinematic Error-Compensation Functions for 

25 Volumetric Error Measurements", Annals of the CIRP, Vol. 34/1/1985, pages 435-438; D.R. 
Hush et al., "An Overview of Neural Networks", Informatica y Automatica, vol 25, 1992; T. 
Moriwaki, C. Zhao, "NN Approach to Identify Thermal Deformation of Machining Center", 
Human Aspects in Computed Integration; and J.C. O'Brien; J.R. Leech, "Can Neural Nets 
Work in Condition Monitoring", Comadem 92 , 1992, pages 88-93). 
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The result provided by such techniques is a computation engine that, given the desired 
model, acquires measures till the model is solved by discarding redundant measures or 
measures that create bad problem conditionings and solution instability. 

This process allows extending the patent applicability to every kind of configuration of 
head 1 and/or head with table. The model can thereby arrive to such complexities as to also 
take into account possible positioning errors of the ball 20 of the gage tool 16 deriving from 
movement errors of the linear axis (those that move the head 1). 

If the error compensation software that resides in the CNC has not the same degree of 
complexity of the measuring software, the same measuring software would perform the 
translation of parameters in its model into the parameters of the compensation software 
module. The requirement of such translation is always obtaining the least degree of residual 
errors in the head 1 . 

EXAMPLE 4: 

As regards the measuring process for the errors that cannot be modeled, all errors not 
included in the geometric model being used are still a reason for positioning errors as regards 
programmed movements of the head 1 axis. 

For such errors it is thereby possible to establish only an empirical link with the 
position of the head 1 axis. The relationship will thereby be univocal. Given a position of the 
two half-bodies A and B, there are three error values DX, DY and DZ in positioning the tool 
tip. The only algebraic relationship is linked to the tool length, such relationship being able to 
be established only by identifying the three values DX, DY and DZ, for a given A and B, with 
two tool length values (therefore the measures for every head 1 position are two, with two 
gage tools 16 with a known and distinct length). For a generic tool length, the error is a linear 
interpolation of the pair of three measured values DX, DY and DZ. 

The measuring process is therefore measuring errors DX, DY and DZ of the ball 20 of 
the gage tool 16 for all affected positions A and B. The measures are carried out first with a 
gage tool 16 and then with the following gage tool 16. The affected positions could be all 
combinations of positions of the two axes from their negative end to their positive end, 
discretizing with a step to be empirically found or with analysis algorithms for error 
frequencies. 



-16- 



The measure of errors that cannot be modeled is described and solved more widely in 
document WO-A-00/003312 of the same Applicant of the present invention. In fact, the 
measure of errors that cannot be modeled is not per se an innovation of the compensation 
method (already described in the above document), but instead in this case the application of 
sensors 14 and support base 11 is relevant in order to automatically perform also this type of 
measures. 

It is necessary to underline that a similar measuring process can be performed by using 
a digital feeler pin (assembled on a spindle nose) and a reference ball (placed in the machine 
working area). 

The parameter input table for such static error compensation software will be composed 



as follows: 
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The error compensation software of the head 1 resides in the CNC. It is also divided 
into a part that performs the head movements compensation using a geometric model and a 
part that is the software that performs the compensation of static, that is "non geometric" (not 
able to be modeled) errors. 

The compensation algorithms are not characterized by any degree of innovation since 
they are already known and they will not be further described. 

As regards the dynamic head test, the above-described system of the invention can be 
used for checking the dynamic behaviour of the controlled axes of the heads 1 and/or the 
movements of these axes combined with those of the linear axis. 

By controlling through the control system 34 the movements to single machine axes, 
the real response on the tool tip can be observed similarly with other control means (for 
example the KGM system manufactured by Heidenhain). 
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5 



10 



Moreover, by controlling the combined movements between linear and rotating axes, 
adapted to obtain that the tool tip (now the ball 20) keeps its position in the space even during 
acceleration transients, the real response curve on the tool tip can be recorded. 

The response curve analysis of the axes compared with the programmed one would 
allow analyzing typical problems like "twitching" and "ripple". 

The only requirement in order to be able to perform this type of test is a pass-band that 
is enough for the signal sampled by the system of the invention. The related pass-band can be 
evaluated in about 1 kHz. 

It has thereby been disclosed that the system of the invention is able, through a single 
measure obtained by the sensors 14 and related to the distance that separates the sensors 14 
themselves from the ball 20, to detect the XYZ coordinates of the center of a tool in a position 
of interest, instead of having, like in the prior art, to perform inaccurate measures, or 
measures that are not located in the XYZ space, or also measures that do not reconstruct the 
tool position in the real position of interest. 

EXAMPLE 5: 

The herein below described solution is the description of a practical arrangement 
resulting from joining together the concepts disclosed in Examples 1, 3 and 4. 
The following components have been used in this example: 

• HMS: measuring system composed of 3 distance sensors that provide an 
incremental digital signal from which their position is directly derived. Such system 
must be placed on the machine tool table that has the head to be measured. 

• Test cylinders: 2 test cylinders with ball end whose length and diameter are known. 
Such cylinders are assembled, through a suitable adapter, on the head spindle. 

• Measuring SW: the software that performs the measuring procedures and computes 
the corrections. 

• Herein below, as an example, reference will be made to a fork-shaped head where 
the first rotating axis is C and rotates around Z while the second rotating axis is A 
and rotates around X when C is 0. 

• Axis strokes will be: 

• axis C from -180 to 180 degrees, 

• axis A from -90 to 90 degrees. 
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The sensors are arranged along the edges of a tetrahedron. Then, seen in plan view, 
they form an angle of 120 degrees one to the other while then they form an angle of 55 degrees 
with respect to the horizon. 

The reason for the tetrahedron arrangement is avoiding collisions between the sensors 
5 and the test cylinder at least in all slants of the positive ball overhanging the device. 

The reason for the 55 degrees is avoiding collisions also with the spindle to which the 
test cylinder is connected. A greater slant would remove further collision problems, thereby 
providing a wider measuring range to the system, while a lower slant tends to better balance 
the measuring resolution in orthogonal directions XYZ. 55 degrees is thereby evaluated as the 
10 best compromise. 

The system measuring position is the one for which the test cylinder ball 
simultaneously touches the three sensors, reaching for each one of them the value for which 
the axis of the three sensors simultaneously passes through the ball center. Consequently in 
such position the tip of every sensor touches the ball in a tangent point, and the plane 
15 perpendicular to the sensor axis is tangent to the ball. To simplify, it is assumed that the 
system is arranged so that the tetrahedron described by the three sensors has its base on the 
plane XY and the test cylinder is directed in a negative direction with respect to Z when A and 
C are 0. 

The system calibration consists in determining: 
20 • the previously described measuring position, and 

• the directing cosines of the three sensors. 

The measuring software, for every head position, has available: 

• the absolute position of rotating axes (C, A) provided by the axis measuring system, 

• the absolute position of linear axes (X, Y, Z) provided by the axis measuring 
25 system, 

• the absolute position of the three sensors (SI, S2, S3). 

In order to determine the measuring position (from now on called "position 0" or 
"sensors position 0"), the software requires that: 

• the ball is previously taken in contact with the sensors, 
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• the position value of the three sensors is given, for which such position is defined 
(this is a construction data related to system construction geometry referred to ball 
diameter). 

After this: 

5 1 .) a movement Z is carried out along a negative direction till at least one of the three 

sensors has reached the center value; 

2. ) from the reached position, a movement X is carried out and, depending on values 
recorded at movement start and arrival, for every sensor the gain is determined with respect to 
movements along machine direction X. Namely, the ratio is determined between value 

10 variation of every sensor and X variation at the movement ends; 

3. ) after having returned to the starting point (the one determined at the end of 
movement 1), an operation similar to step 2 is carried out to determine, for every sensor, the 
gain with respect to movements along direction Y, then the movement is carried out along 
direction Y; 

15 4.) a series of movements XY are carried out, adapted to search for the values SI, S2, 

S3 that better approximate position 0. Since probably position Z is not the correct one, the 
searched values are those for which this position is defined apart from a scale factor. Then, 
depending on previously computed gains, step by step there are determined movements XY, 
which every sensor would require to assume value 0. It must be taken into account that the 

20 search was successful if the values are determined apart from a certain tolerance. If the search 
does not converge, it is stopped and restarted from step 1. This occurs when the starting 
position was for the step 1 to load the sensors in an unequal way so that steps 2 and 3 
determined very different gains with respect to the ones in point 0. The first iteration in step 4 
anyway has allowed to adequately approach such condition. Then, the new execution of steps 

25 1,2 and 3 provides stronger values that will take step 4 to converge; 

5.) if step 4 has reached the convergence, a movement Z is again carried out since step 
4 has determined position 0 apart from a scale factor for SI, S2 and S3! 

In order to determined directing cosines for the three sensors, starting from the 
previously determined position 0, a series of movements XYZ are carried out, adapted to 

30 determine, for every sensor, the relationship (gain) between sensor measure and movement of 
one of linear axes. This is determined with search algorithms that are based on the fact that, 
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for position 0, where the sensors axis passes by the ball center and the plane perpendicular 
thereto is tangent to the ball, the following is valid: 

• describing a curve that has as abscissa the angle, in plane XY, formed by a 
movement XY that starts from position 0 and as ordinate the elongation value of a 
sensor for such movement, such curve has a single maximum coinciding with the 
movement XY whose angle coincides with the projection in XY of the sensor axis; 

• having determined through the previous step the direction in XY of every sensor, a 
similar principle is valid for the curve described by movements XYZ for which the 
straight line along which the movement is carried out lies in the plane perpendicular 
to XY and oriented according to the angle XY determined in the previous step. The 
curve has, as abscissa, the angle formed by movements XYZ with respect to plane 
XY and as ordinate the elongation value of a sensor. The curve maximum points 
out in this case the movements XYZ in which the performed distance coincides with 
the elongation that the examined sensor had. The ratio between sensor elongation 
and movements X, Y and Z for this last position are the searched gains. 

The measure principle is based on the fact of being able to drive the ball to position 0. 
For every sensor the following are known: 

• directing cosines Kx#, Ky#, Kz#, where # represents the current sensor (1, 2, 3), 

• the sensor value in position 0 (S0#). 

For a generic ball position, described by XYZ and SI, S2 and S3, the following are 
determined: 

Movement X required by sensor # = MX# = (S0# - S#) x Kx# 

Movement Y required by sensor # = MY# = (S0# - S#) x Ky# 

Movement Z required by sensor # = MZ# = (S0# - S#) x Kz# 

Then the required movement will be: 

Axis movement X = MX1 + MX2 + MX3 

Axis movement Y = MY1 + MY2 4- MY3 

Axis movement Z = MZ1 + MZ2 + MZ3 

These movements are then carried out by a numeric control system that continuously 
reads the values S0# - S# and inserts the movements to be carried out, performing suitable 
acceleration and deceleration logics. The algorithm then operates as if the three sensors are 
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three springs that, connected to the ball, tend to take back the ball to a balanced position, 
defined apart from a certain tolerance, that in principle coincides with position 0 but must not 
compulsorily do so. In fact, in order to best calibrate the search algorithm for the balanced 
position, the values S0# are placed equal to the value that the sensors assume in their balanced 
5 position as soon as it is possible to reach it for the first time. 

Then, in practice, the calibration procedure, that determines point 0, for which 
S#=S0#, and directing cosines Kxyz#, has not the purpose of exactly determining the HMS 
arrangement but the purpose of providing to the previously described algorithm, data that 
allow always determining a balanced position for the ball. In practice this balanced position 
10 can represent, in the Cartesian space, an always identical position of the ball, that modifies the 
sensors contact according to a possible ball rotation, and to approximations and computation 
tolerances of the algorithm, that, being a system for successive and continuous approximations, 
will tend to minimize them. 

The major advantages of this procedure are: 
15 • the search for the balanced position is minimally affected by the accuracy with 

which S0# and Kxyz# have been determined, and thereby the calibration algorithms 
are particularly simple and quick, 

• no computation complexity is added if the S0# of the three sensors are different as 
well as if the sensors are not exactly oriented along the sides of the ideal 

20 tetrahedron. This allows realizing the HMS without particular restrictions on 

dimensional tolerances both of components and of assembling, and allows placing 
the HMS on a machine tool without particular care for its orientation, 

• if during the measures, thermal effects distort the HMS morphology, the influence 
on Kxyz# values has no reflexes on search stability for the balance point and 

25 therefore for the measures that have to be carried out. Therefore it is useless to 

carry out again the calibration and the system is particularly stable, 

• the search for the balance position is quick and therefore the measure of an 
individual point is particularly efficient. 

These advantages are difficult to be obtained with any other system whose purpose is 
30 identifying the relationship that directly links, for a given triad of sensor values, the 
corresponding position value XYZ. 
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Then, it is now possible to deem that it is safe to ask the system to take back the ball to 
such a position that this position will always be the same, if referred to an inertial Cartesian 
system, apart from: 

• ball sphericity error (the ball can be easily obtained with sphericity tolerances lower 
than 2 micrometers), 

• position accuracy for machine axes XYZ, 

• computation accuracy (practically neglectable), 

where the greatest, and practically the only, influence is given by machine positioning 
accuracy. 

It will be possible to state that the system "centers the ball" when it carries out such 
operation. 

A correct calibration is easily checked by performing a series of positionings XYZ and 
verifying that, starting from such positioning, the reached balance point is always the same 
point XYZ apart from a certain tolerance. 

As regards the measure of a point, it must be remembered that the measuring software 
has available, for each head position: 

• the absolute position of rotating axes (C, A) provided by the axis measuring system, 

• the absolute position of linear axes (X, Y, Z) provided by the axis measuring 
system, 

• the absolute position of the three sensors (SI, S2, S3). 

In practice the system has access to different XYZ values, of which: 

• the absolute position of machine axes directly provided by the axis measuring 
system (XA, YA, ZA, AA, CA) (in MXA1 CNC manufactured by Fidia S.p.A., 
Italy), 

• the absolute position of machine axes provided by the control system after having 
applied a correction to the one provided by the axis measuring system (X, Y, Z, A, 
C); this transformation, performed by the control, is adapted to correct the values 
provided by the measuring systems in order to approach the XYZ system to an 
inertial Cartesian reference system (in Fidia MXA1 CNC), 

• the relative ball center position (XL, YL, ZL); depending on XYZ AC values, on 
geometric ball shaping data and on test cylinder dimensions, the numeric control 
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provides the supposed ball center position (in Fidia MXA1 CNC, POSITION with 
RTCP ON). 

Error measures are then carried out according to the following procedure: 

• for a given head position 1 (Al, CI), the ball is centered, and absolute and relative 
5 XYZ values are recorded (XI, Yl, Zl, XL1, YL1, ZL1), 

• the head is moved to position 2 (A2, C2), the ball is centered, and absolute and 
relative XYZ values are recorded (X2, Y2, Z2, XL2, YL2, ZL2), 

• variations of the 2 positions are computed: 
Dx = XI - X2 

10 Dy = Yl - Y2 

Dz = Zl - Z2 
DLx = XL1 - XL2 
DLy = YL1 - YL2 
DLz = ZL1 - ZL2 

15 Then, for the movement from position 1 to position 2, Dxyz represents, in the 

Cartesian reference system, the movement that the head must perform in order not to make 
ball positioning errors, while DLxyz represent the error that is made by such movement. 
DLxyz are then the error measure provided by the sensors as if they were oriented according 
to a Cartesian triad. 

20 Movement from position 1 to 2 is carried out by the control taking care that the ball 

does not abandon the contact with sensors (RTCP function), otherwise it would not be possible 
to center again the ball. 

As regards measuring procedures, it has been pointed out that, given two programmed 
head positions, the system determines the ball center positioning error and reference system 
25 positions XYZ if such errors are null. 

Typically, the measuring functions perform two types of measures: 

1 . ) simple measures that determine Dxyzac and DLxyzac from two or more measuring 
positions, 

2. ) complex measures in which a rotating axis is not moved while the other one 
30 performs a series of positionings from a starting position to an end position by performing 

prefixed steps. By carrying out such measure, the ball, as said, remains ideally unmoving, and 
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the linear axes perform a circle in space. This type of measure will be called the "measure of a 
circle". The measure of a circle will have the following notation: 

() [Moved axis: C/A] ; [starting position in degrees - START] ; [end position in 
degrees - END] ; [pitch in degrees - STEP] 
5 ® [Fixed axis: A/C] ; [Fixed axis position in degrees] 

From these measures, lists of values will be recorded that are Dxyzac and DLxyzac of 
every measure point. 

Known computation algorithms are then able to determine, from measures of a circle, 
geometric entities such as: 
10 • center and radius of the circle described by Dxyzac, 

• plane passing through Dxyzac, 

• cube curve passing through Dxyzac for a given quadratic error tolerance. 
From many measures of circles and from related obtained geometric entities, it is 

thereby possible to determine values and other geometric entities, such as: 
15 • distance between planes, 

• intersections of planes and of planes with cubic curves, 

• cones and/or cylinders passing through many circles, 

• straight lines tangent to many circles, 

• straight lines passing through the center of many circles, etc. 

20 The angular head positioning error is corrected by inserting corrections that allow the 

control to determine A and C given AA and CA. 

The computation of such corrections is based on the principle that, after having 
performed a circle for A or C, the angle, that is determined by joining the center of a circle 
computed with individual measuring points, represents, apart from a constant, the remaining 

25 correction to be inserted to pass from A A and/or CA to A and/or C. 

The above mentioned constant depends on the influence of distances between head 
rotation axes (head misalignments) and the fact that the rotation axes zero position has still not 
been correctly defined. This constant is removed if these measures are followed for two mirror 
positions of the axis that is not moved. By averaging the values obtained by the 2 circle 

30 measures, the constants mutually cancel because their signs are opposite. 
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The two arrangements further allow taking into account and averaging also angular 
errors inserted by the different weight distributions that the head could assume in the two 
arrangements. 



As regards the measuring procedure for the angular positioning error on axis C, such 
measure allows performing two series of circles for mirror positions of axis A. 



Circle 1 


Circle 2 


() C ; -180; 180; [STEP] 
® A ; +90 


O C ; -180; 180; [STEP] 
® A ; -90 



Axis A is slanted at its maximum stroke to obtain the highest possible resolution in 
computations, namely the circle with maximum radius. If the head has such misalignments that 
the two circles have different radius, the two positions of A must be chosen also in order to 
make the circle radius equal in order to keep the computation resolution constant. 



Every circle must be performed many times in order to make averages on measures 
both along a positive direction and along a negative direction, also in order to determine the 
inversion error of the examined rotating axis. 

The difference between programmed position for axis C and computed angle for the 
straight line joining the circle center and a single measured point XYZ, is the remaining 
correction to be added to pass from CA to C. 

Data computed between the two circles are averaged in order to: 

• remove the effects of constants, 

• take into account angular errors inserted by the weight that could be distributed 
differently with A at 90 and -90 degrees. 

The final results are anyway translated in order to minimize inserted correction values. 
Such translation will then be taken into account and corrected during the zero point measure 
for axis C. 

In order to perform this measure, no other parameter must be known. 

As regards the measuring procedure of the angular positioning error on axis A, such 
measure consists in performing two series of circles in order to take into account errors both in 
XZ and in YZ planes. 
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Circle 1 


Circle 2 


O A ; -90; 90; [STEP] 
® C ; 0 


OA; -90; 90; [STEP] 
® C ; 90 



Every circle can be performed many times in order to average the measures both along 
a positive direction and along a negative direction, in order to also determine the inversion 
error of the examined rotating axis. 



The difference between programmed position for axis A and angle computed for the 
straight line joining the circle center and a single measured point XYZ, is the remaining 
correction to be added to pass from A A to A. 

Data computed between the two circles are averaged in order to: 

• remove the effect of constants, 

• take into account angular errors inserted by the different weight distribution with C 
ad 90 and 0 degrees. 

The final results are anyway translated in order to minimize the inserted correction 
values. Such translation will then be taken into account and corrected during the zero point 
measure for axis A. 

In order to carry out this measure, no other parameter must be known. 

On all measured data, the following computations are made: 

• concentricity of computed circles, 

• mean, standard deviation and band of data differences, 

• repeatability on the same position, 

• in order to determine the efficiency of collected data. 

The rotating axes zero points for a head AC are those for which: 

• with C = 0 axis A is moved on plane YZ 

• with A = 0 the test cylinder is perpendicular to plane XY 
Such values are commonly called "zero preset" of a rotating axis. 

The zero preset measure of C consists in performing four series of circles in order to 
take into account both XZ and YZ planes along many rotation directions of axis A. 
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Circle 1 


Circle 2 


OA; -90; 90; [STEP] 
® C ; 0 


OA; -90; 90; [STEP] 
® C ; 90 


Circle 3 


Circle 4 


() A ; -90; 90; [STEP] 
® C ; 180 


OA; -90; 90; [STEP] 
® C ; -90 



From circles 1 and 3 a plane is determined and therefore the angle that it forms with 
plane YZ. 

From circles 2 and 4 a plan is determined and therefore the angle that it forms with 
5 plane XZ. 

The best trade-off for C preset is then searched with averages in order to align the 
above mentioned planes and then find the best head quadrature. 

To carry out this measure, no other parameter must be known. 

The zero preset measure for A consists in performing four series of C circles, all with 
10 axis A at 0. Every circle is performed with a different test cylinder arrangement. 



Circle 1 


Circle 2 


<) C ; -180; 180; [STEP] 
® A ; 0 

Cylinder with length LI in position 1 


<> C ; -180; 180; [STEP] 
® A ; 0 

Cylinder with length LI rotated by 180 degrees 
with respect to 1 


Circle 3 


Circle 4 


() C ; -180; 180; [STEP] 
® A ; 0 

Cylinder with length L2 in position 1 


O C ; -180; 180; [STEP] 
® A ; 0 

Cylinder with length L2 rotated by 180 degrees 
with respect to 1 



The average between circles 1 and 2 and then 3 and 4 allows removing possible errors 
inserted by non-rectilineity of the test cylinder. 

The two found circles describe a cone, also knowing position Z of the two centers, or 
the difference between LI and L2. 
15 The cone angle will be the value to be corrected for axis A preset. 
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The greater L1-L2 is, the greater the computation resolution will be. 
In order to carry out this measure, no other parameter must be known. 
As regards the head geometry measure, for a head AC the geometric values are as 
follows: 

Misalignment name Description 

1) DXBT With C=0 and A=0, distance along direction X between C 

rotating axis and test cylinder axis. 

2) DZBA With C=0 and A=0, distance along direction Y between C 

rotating axis and A rotating axis. 

3) DZAT With C=0 and A=0, distance along direction Y between A 

rotating axis and test cylinder axis. 

4) RTCPLKS With C = 0 and A = 0, distance along direction Y between A 

rotating axis and test cylinder connection point (point from 
which the cylinder length is then determined). 

5) RTCPANG With C=0 and A=0, angle in plane XZ between C rotating 

axis and A rotating axis. 

6) With C=0 and A=0, angle in plane YZ between C rotating 

axis and A rotating axis. 

7) Angle in plane XZ between C rotating axis and axis Z. 

8) Angle in plane YZ between C rotating axis and axis Z. 
5 6, 7, and 8 cannot be corrected with Fidia controls. 

The DXBT measure consists in performing four series of circles in order to take into 
account both XZ and YZ planes along many rotation directions of axis A. 
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Circle 1 


Circle 1 


OA; -90; 90; [STEP] 

® C ; 0 


O A ; -90; 90; [STEP] 
® C ; 90 


Circle 3 


Circle 4 


OA; -90; 90; [STEP] 
® C ; 180 


OA; -90; 90; [STEP] 
® C ; -90 



The distance between planes for circle 1 and circle 3, as well as the distance between 



planes for circle 2 and circle 4 are twice the DXBT. 

DXBT can further be determined by searching the value that reduces at a maximum the 
5 difference between: 

• XL of circle 1 and XL of circle 2, 

• YL of circle 3 and YL of circle 4. 

The found value will be the remaining correction to be given to DXBT. 
Usually, when possible like in this case, it is more efficient to determine the value of a 
10 parameter depending on effects that it produces on ball positioning errors (namely the various 
XYZL) rather than starting from geometric head data (namely from XYZ). 

In order to carry out this measure, the axis C preset must be known, after having 
previously measured it. 

The DZBA measure performs a circle of C with A at 0. 
15 () C ; -180; 180; [STEP] 
® A ; 0 

The circle radius will be: 
RADIUS = RADQ( DXBT"2 + (DZBA + DZAT)"2 ). 

Being then DXBT known, the sum of DZBA and DZAT is determined. 
20 As first approximation, DZAT will be considered as null and DZBA only is corrected. 

Afterwards, the measure to compute DZAT will correct the performed approximation. 
In order to carry out this measure, the following must be known, after having 
previously been measured: 
© DXBT, 
25 • axis A preset. 
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The DZAT measure consists in performing four series of circles in order to take into 
account both XZ and YZ planes along many rotation directions of axis A. 



Circle 1 


Circle 2 


O A ; -90; 90; [STEP] 
® C ; 0 


O A ; -90; 90; [STEP] 
® C ; 90 


Circle 3 


Circle 4 


() A ; -90; 90; [STEP] 
® C ; 180 


() A ; -90; 90; [STEP] 
® C ; -90 



At this time, DZAT is null and therefore the control does not perform therefore any 
compensation of errors YL of circles 1 and 3, while errors XL of circles 2 and 4 will be its 



5 image. 

Such errors describe a circle since the compensation due by DZAT would be maximum 
for A=0 and null for A = 90 and A = -90. 

The radius that better approximates this circle is the best compromise for DZAT. The 
same value must then be subtracted from DZBA following the approximations inserted by its 
10 computation. 

In order to carry out this measure, the following must be known, after having 
previously measured them: 

• DXBT, 

• DZBA, 

15 • axis A preset. 



The RTCPLKS is determined as average of radius of the following circles: 



Circle 1 


Circle 2 


() A ; -90; 90; [STEP] 
® C ; 0 


C) A ; -90; 90; [STEP] 
® C ; 90 


Circle 3 




() C ; -180; 180; [STEP] 
® A ; 90 





To carry out this measure, everything must be known, after having previously 
measured it. 



As regards the RTCPANG measure, the projection in XZ of the circle plane vector: 
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OA; -90; 90; [STEP] 
® C ; 0 

is the RTCPANG value. 

Regarding the test cylinder, its length (measured from the point where the RTCPLKS 
5 measure ends till the tip of the cylinder itself) and the ball diameter must be known. 

Once having assembled the cylinder in the spindle, due to its shape or due to 
assembling, the ball center could have deviations from the spindle axis. 

Such deviations generate a direct error in measuring the DXBT and DZAT parameters. 
Once having known these deviations, it is possible to carry out the above described 
10 measures and at the end simply subtract these values from DXBT and DZAT. 



• half of the difference between values of the two measures are the searched values. 

Since the manufacture of the test cylinder and the adapter to assemble it on the spindle 
tends to minimize at a maximum deviations of the ball center with respect to spindle axis, 
measures along directions X and Y will be small values, while the measure along direction Z 
20 will be their consequence and by at least two order of magnitude lower. The error along 
direction Z will then be neglectable. 

Since the previous measure is affected by the accuracy with which the 180-degree 
rotation is performed, another method to determine the same values is: 



15 



In order to measure these values, the following procedure can simply be carried out: 

• assemble the test cylinder in the spindle 

• carry out the ball centering and record XL, YL and ZL, 

• rotate the cylinder by 180 degrees in the spindle, 

• carry out the ball centering again and record XL, YL and ZL, 



assemble the test cylinder in the spindle, 



25 



carry out a ball centering and record XL, YL and ZL, 



repeat the same operation by rotating the cylinder at 90-degree pitches and 
recording XL, YL and ZL, 

by carrying out the rotation and measuring operations four times, the angular 
positions in which it is assumed to have rotated the cylinder will be: 



30 



• starting position 0, 

• position 90, 
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• position 180, 

• position 270, 

• return to starting position. 

In this case it is assumed that one has been able to obtain angular positions approximate 
5 to ± 5 degrees. In this case the problem is determining: 

• the circle that passes through the list of measured points XL and YL. Since the 
points are five, the circle and its radius can be obtained with the minimum square 
method, 

• determine the final angle at which the cylinder is placed. Being circle center and 
10 radius known in the same coordinate system of XL and YL, such value can be 

determined algebraically and the circle radius can be correctly decomposed in test 
cylinder misalignment values. 
It is clear then that the cylinder must not be moved any more once having located its 
misalignment values unless one has again to carry out the calibration. 
15 As regards the dependency on measures, X points out the values that affect the 

measures and that therefore must be previously be measured. 

V points out the values that affect the measures but that are kept into account by the 
measuring procedure itself and then are removed during computations. 



Dependencies refer to the previously shown measuring procedures. 



Measure of: 


DXBT 


DZAT 


DZBA 


RTCPLKS 


PRESET A 


PRESET C 


T.c.c.* 


Angular position C 












V 




Angular position A 














V 


PRESET C 
















PRESET A 














V 


DXBT 












X 


X(*) 


DZBA 


X 








X 






DZAT 


X 




X 




X 




X(*) 


RTCPLKS 


X 


X 


X 




X 




X 


RTCPANG 












X 





20 T.c.c.: Test Cylinder Coaxiality 
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(*) : Like for other measures, the X-type influence presumes the knowledge of the 
influence value in order to correctly determine the measured parameters. Otherwise, for these 
cases, it is noted that the dependency must be solved a posteriori. For example, while for 
DXBT the PRESET C must be known, measured and corrected in order to correctly perform 
5 the DXBT computation algorithm, instead in the DZAT case, the non-coaxiality cylinder value 
has no influence on computation algorithms since the measure result will in practice be a value 
that is equal to DZAT + cylinder non-coaxiality. Then the correct DZAT value will be: 
measuring and computation algorithm result - cylinder non-coaxiality. The same holds for 
DXBT, always referred to test cylinder non-coaxiality. 
10 The system is also able to carry out an approximate error check for axes XYZ and also 

AC that will not be able to correct since they cannot be included in head corrections. 

As regards the rotation planarity check for axis C and the rectilineity check for linear 
axes XY, the following circles are performed: 



Circle 1 


Circle 2 


<) C ; -180; 180; [STEP] 
® A ; 90 


() C ; -180; 180; [STEP] 
® A ; -90 



Taking into account that position A-90C180 corresponds to position A90C0, as well as 
15 position A-90C-90 corresponds to position A90C90 and so on, data are ordered and the 
following is computed: 

• the circle that passes from even points at half of differences among values of X, Y 
and Z. The vectors of such circle plane represent the rotation planarity errors of 
axis C. 

20 • the curves for even points at half of the sum of values of X, Y and Z. Such curves 

represent the rectilineity errors of axes X and Y along direction Z. 
As regards the orthogonality check for axes XYZ, the following circles are performed: 



Circle 1 


Circle 2 


<) C ; -180; 180; [STEP] 
® A ; 90 


O C ; -180; 180; [STEP] 
® A ; 0 



The difference in X and Y of the circle centers represents the orthogonality error of 
axis Z with respect to X and Y. 
25 As regards the axis scale check, the following circle is performed: 
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O C ; -180; 180; [STEP] 
® A ; 90 

For such series of data, the following are determined: 
• the radius of the circle that best approximates them, 
5 • the axes of the ellipse that best approximates them. 

The difference between circle radius and ellipse axes is a scale error in data provided 
by XY measuring systems. 

As regards the rotation check for axis C, the following circle is performed: 
(> C ; -180; 180; [STEP] 
10 ® A ; 90 

The offset of the curve of Z values from a simple sinusoid period (that would be the 
effect of C non-planarity) is the rotation axis swinging for axis C. 

As regards the rotation check for axis A and the rectilineity check for linear axes XYZ, 





the following circles are performed: 






Circle 1 


Circle 2 




C) A ; -90; 90; [STEP] 
® C ; 0 


() A ; -90; 90; [STEP] 
® C ; 180 


15 


Taking into account that position A-90C1! 


30 corresponds to position A90C0, as well as 



position A90C180 corresponds to position A-90C0 and so on, data are ordered and the 
following is computed: 

• from half of the differences of values of X, Y and Z the rotation errors for axis A, 

• at half of the sums of values of X, Y and Z the rectilineity errors of axis Y along 
20 direction X. 

By instead performing the circles: 



Circle 1 


Circle 2 


OA; -90; 90; [STEP] 
® C ; 90 


O A ; -90; 90; [STEP] 
® C ; -90 



the same values are obtained for axis A and furthermore the rectilineity errors of axis X 
along direction Y. 

By finally composing the data of the four circles, it is further possible to extract from 
25 rotation errors for axis A, the rectilineity influence of axis Z. 
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The measure of remaining errors collects data in XL, YL and ZL for all positions of A 
and C that describe their whole strokes for a predefined pitch. In practice the XYZ value 
would describe a ball portion. 

This is carried out with two test cylinders with different length and by rotating them by 
5 180 degrees to remove, performing their averages, the influences of its lack of axiality. Then 
the sets of performed measures are four. 



The collected data are used to fill-in the following table: 



Position A 


Position B 


Bar length 


DX 


DY 


DZ 




[negative end of 
stroke] 


S&ig&lfl 


# 


ft 


tt 


oooooo 


[negative end of 
stroke] 


ILsBg&fl 


tt 


ft 


ft 




[negative end of 
stroke] in 




# 


it 


ft 


[negative end of 
stroke] 




ILas!g& A 


# 


it 


ft 







ILasgiii a 




it 


ft 


[positive end of 
stroke] 






# 


it 


ft 


[negative end of 
stroke] 


[positive end of 
stroke] 




if 


it 


tt 




[positive end of 
stroke] 




# 


ft 


tt 


[positive end of 
stroke] 


[positive end of 
stroke] 


ILsag& fl 


ft 


tt 


ft 


[negative end of 
stroke] 


[negative end of 
stroke] 




# 


ft 


ft 




[negative end of 
stroke] 




ft 


ft 


ft 


[positive end of 


[negative end of 




ft 


tt 


tt 
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Position A 


Position B 


Bar length 


DX 


DY 


DZ 


stroke] 


stroke] 




3 








stroke] 






U 
ft 


u 

if 


U 

ff 








% 

% . 


U 
ft 


U 
It 


u 

ft 


[positive end of 






tt 


# 


tt 


[negative end of 

O LI v/rv^ J 


[positive end of 

O LI UA^J 




% 


tt 


tt 


tt 




[positive end of 
stroke] 


Bge'n'gtB 


% 


tt 


tt 


tt 


[positive end of 
stroke] 


[positive end of 
stroke] 


RgengtE 


% 




tt 


tt 



The numeric control searches and interpolates in such table the values that allow it to 
compute the correction along direction XYZ for a given position of A and C and a given tool 
length. 



Although the present invention has been discussed in considerable detail with reference 
to certain preferred embodiments, other embodiments are possible. Therefore, the scope of 
the appended claims should not be limited to the description of preferred embodiments 
contained in this disclosure. All references cited herein are incorporated by reference to their 
entirety. 
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